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Abstract 

Dense random packings of hard particles are useful models of granular media and are closely 
related to the structure of nonequilibrium low-temperature amorphous phases of matter. Most work 
has been done for random jammed packings of spheres, and it is only recently that corresponding 
packings of nonspherical particles (e.g., ellipsoids) have received attention. Here we report a 
study of the maximally random jammed (MRJ) packings of binary superdisks and monodispersed 
superballs whose shapes are defined by \x\\ 2p + ■ ■ ■ + \xd\ 2p < 1 with d = 2 and 3, respectively, 
where p is the deformation parameter with values in the interval (0, oo). As p increases from zero, 
one can get a family of both concave (0 < p < 0.5) and convex (p > 0.5) particles with square 
symmetry (d = 2), or octahedral and cubic symmetry (d = 3). In particular, for p = 1 the particle 
is a perfect sphere (circular disk) and for p — > oo the particle is a perfect cube (square). We find 
that the MRJ densities of such packings increase dramatically and nonanalytically as one moves 
away from the circular-disk and sphere point {p = 1). Moreover, the disordered packings are 
hypostatic, i.e., the average number of contacting neighbors is less than twice the total number of 
degrees of freedom per particle, and the packings are mechanically stable. As a result, the local 
arrangements of particles are necessarily nontrivially correlated to achieve jamming. We term such 
correlated structures "nongeneric" . The degree of "nongenericity" of the packings is quantitatively 
characterized by determining the fraction of local coordination structures in which the central 
particles have fewer contacting neighbors than average. We also show that such seemingly "special" 
packing configurations are counterintuitively not rare. As the anisotropy of the particles increases, 
the fraction of rattlers decreases while the minimal orientational order as measured by the cubatic 
order metric increases. These novel characteristics result from the unique rotational symmetry 
breaking manner of the particles, which also makes the superdisk and superball packings distinctly 
different from other known nonspherical hard-particle packings. 

PACS numbers: 61.50.Ah, 05.20.Jj 
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I. INTRODUCTION 



Particle packing problems, such as how to fill a volume with given solid objects as densely 
as possible, are among the most ancient and persistent problems in science and mathematics. 
A packing is a large collection of non-overlapping solid objects (particles) in <i-dimensional 
Euclidean space M. d . The packing density <fi is defined as the fraction of space M. d covered by 
the particles. Dense ordered and random packings of nonoverlapping (hard) particles have 
been employed to understand the equilibrium and non-equilibrium structure of a variety 
many-particle systems, including crystals, glasses, heterogeneous materials and granular 
media jl-4]. Packing problems in dimensions higher than three attract current interest for 
retrieving stored data transmitted through a noisy channel 5|48]. 

The packings of congruent hard spheres in M 3 have been intensively studied since despite 
the simplicity they exhibit rich packing characteristics. It is only recently that the densest 
packings with (p max = 7r/vT8 ~ 0.74, realized by the face-centered cubic lattice and its 
stacking variants, have been proved 9|. In addition, three-dimensional random packings 
can be prepared both experimentally and numerically with a relatively robust density <ft ~ 



0.64 



10 



111 ]. The term random close packing (RCP) 12], widely used to designate the 



"random" packing with the highest achievable density, is ill-defined since random packings 
can be obtained as the system becomes more ordered and a definition of randomness has 
been lacking. A more recent concept that has been suggested to replace RCP is that of 
the maximally random jammed (MRJ) state ll|, corresponding to the most disordered 
among all jammed (mechanically stable) packings. A jammed packing is one in which the 
particle positions and orientations are fixed by the impenetrability constraints and boundary 
conditions [l^j]. It has been established that the MRJ state for spheres in M 3 has a density 



of ~ 0.637, as obtained by a variety of different order metrics [13| 



This density value 



is consistent with what has traditionally asscoiated with RCP in three dimensions. 

It has been argued in the granular materials literature that large disordered jammed 
(MRJ) packings of hard frictionless spheres are isostatic flol lit] ], meaning that the total 
number of interparticle contacts (constraints) equals the total number of degrees of freedom 
of system and that all of the constraints are (linearly) independent. This implies that 
the average number of contacts per particle Z is equal to twice the number of degrees of 
freedom per particle / (i.e., Z = 2f), in the limit as the number of particles gets large. This 
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prediction has been verified computationally with very high accuracy [17j, |18] . On the other 
hand, a packing is hypostatic if it is mechanically stable (i.e., jammed) while the number 
of constraints is smaller than the number of degrees of freedom. For large packings, this is 
equivalent to the i neq uality Z < 2f. It has been shown that a jammed sphere packing can 
be not hypostatic [17 ]. 

It is also of great practical and fundamental interest to understand the organizing prin- 



ciples of dense packings of nonspherical particles |19h28]. The effect of asphericity is an 
important feature to include on the way to characterizing more completely real dense granu- 
lar media as well as low-temperature states of matter. Another important application relates 
to supramolecular chemistry 29] of organic compounds whose molecular constituents can 
possess many different types of group symmetries 30]. Such systems can be approximated 
by nonspherical hard particles with the same group symmetries. 
Recently, MRJ packings of three-dimensional ellipsoids 
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32] have been studied. In 



particular, it was found that the density <fi and the average coordination number Z (the 
average number of touching neighbors per particle) increase rapidly, in a cusp-like manner, 
as asphericity is introduced from the sphere point. The density 4> reaches a maximum at a 
critical aspect ratio a* [33] and then begins to decrease; while Z is increasing monotonically 
until it attains the plateau value for all a beyond a*. In addition, Z is always smaller than 
twice the number of degrees of freedom per particle / with its plateau value slightly below 
2f (for an ellipse / = 3 and for an ellipsoid / = 6). In other words, the packings are 
hypostatic 32J. The characteristics of MRJ ellipsoid packings are distinctly different from 
their densest crystalline counterpart 21], in which <fr increases smoothly as one moves away 
from the sphere point, and reaches a plateau value of 0.7707... for a > a/3 (oblate spheroids) 
and a < 1 / \/3 (prolate spheroids) . 

In Refs. [25] and [26], we studied dense and maximally dense packings of superballs, 
a family of nonspherical particles with versatile shapes. In particular, a <i-dimensional 
superball is a centrally symmetric body in M. d occupying the region 



\X!\ 2P + \X 2 \ 2P + ■ ■ ■ + \x d \ 2p < 1, (1) 

where Xi (i — 1, . . . , d) are Cartesian coordinates and p > is the deformation parameter, 
which indicates to what extent the particle shape has deformed from that of a <i-dimensional 
sphere (p = 1). Henceforth, the terms superdisk and superball will be our designations for 
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FIG. 1: (color online). Superdisks with different values of the deformation parameter p. 
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FIG. 2: (color online). Superballs with different values of the deformation parameter p. 



the two-dimensional (d = 2) and three-dimensional (d = 3) cases, respectively. A superdisk 
possesses square symmetry, as p moves away from unity, two families of superdisks can 
be obtained, with the symmetry axes rotated 45 degrees with respect to each other; when 
p < 0.5, the superdisk is concave (see Fig. [T]). A superball can possess two types of shape 
anisotropy: cube-like shapes (for p > 1) and octahedron-like shapes (for < p < 1) with a 
shape change from convexity to concavity as p passes downward through 0.5 (see Fig. [2]). 

Optimal packings of congruent superdisks and superballs apparently are realized by cer- 
tain Bravais lattices possessing symmetries consistent with those of the particles 25|, |26|. |28|. 



Even these crystalline packings exhibit rich characteristics that are distinctly different from 
other known packings of nonspherical particles. For example, we found that the maximal 
density <p m ax as a function of p at p = 1 (the sphere or circular-disk point) is nonanalytic 
and increases dramatically as p moves away from unity. In addition, we have discovered 
two-fold degenerate maximal density states for square-like superdisks, and both cube-like 
and octahedron-like superballs. 

In this paper, we generate both packings of binary superdisks in R 2 and monodisperse 



superballs in M 3 that represent the maximally random jammec 
cles, using a novel event-driven molecular dynamics algorithm 



(MRJ) state of these parti- 
and investigate their 
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characteristics. For both superdisks and superballs, we find that the corresponding density 
4> and the average contact number Z increase rapidly, in a cusp-like manner, as the particles 
deviate from perfect circular disks and spheres, respectively. In particular, we find that the 
MRJ packing density <j) increases monotonically as p moves away from unity, and shows no 
signs of a plateau even for large p values. This is to be contrasted with the case of ellipsoids 
for which the packing density reaches a maximum as the aspect ratio ratio increases from 
its sphere-point value and then begins to decrease as the aspect ratio grows beyond that 
associated with the density-maximum value. 

Moreover, we find that Z for superdisk and superball packings reach its associated plateau 
value at relatively small asphericity deviations (i.e., \p — 1|) and the packings remain hypo- 
static for all values of p examined. By "hypostatic", we mean that the Z is smaller than 
twice the number of degrees of freedom per particle, compared to random ellipsoid packings 
where the plateau value of Z is only slightly below 2f. Therefore, to achieve jamming, 
the local particle arrangements are necessarily correlated in a nontrivial way. We call such 
correlated structures "nongeneric" 36(. We quantify the degree of "nongenericity" of the 
packings by determining the fraction of local coordination configurations in which the cen- 
tral particles have fewer contacting neighbors than average Z. We also show that such 
"nongeneric" configurations are not rare, which is a rather counterintuitive conclusion. In 
addition, we find that although the rapid increase of density is unrelated to a ny observable 
translational order, the orientational order (e.g., the cubatic order parameter 37|) increases 
as p moves away from unity. These packing characteristics, which are distinctly different 
from that of the MRJ packings of ellipsoids, are due to the unique way in which rotational 
symmetry is broken in superdisk and superball packings. 

The rest of the paper is organized as follows: In Sec. II, we briefly describe the simulation 
techniques and the obtained packings. In Sec. Ill, we provide a detailed analysis of the novel 
packing characteristics. In Sec. IV, we make concluding remarks. 



II. MAXIMALLY RANDOM JAMMED PACKINGS VIA COMPUTER SIMULA- 
TION 



We use an event-driven molecu 



Donev, Torquato and Stillinger 
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ar dynamics packing algorithm recently developed by 



35| (henceforth, referred to as the DTS algorithm) to 
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generate MRJ packings of convex superdisks in two dimensions and superballs in three 
dimensions. The DTS algorithm generalizes the Lubachevsky-Stillinger (LS) sphere-packing 
algorithm 38j to the case of other centrally symmetric convex bodies (e.g., ellipsoids and 
superballs). Initially, small particles are randomly distributed and randomly oriented in 
the simulation box (fundamental cell) with periodic boundary conditions and without any 
overlap. The particles are then given translational and rotational velocities randomly and 
their motion followed as they collide elastically and also expand uniformly with an expansion 
rate 7, while the fundamental cell deforms to better accommodate the configuration. After 
some time, a jammed state with a diverging collision rate is reached and the density reaches a 
local maximum value. To generate random jammed packings, initially large 7 are employed 
to prevent the system following the equilibrium branch of the phase diagram that leads to 
crystallization. Near the jamming point, sufficiently small expansion rate is necessary for 
the particles to establish contacting neighbor networks and to form a truly jammed packing. 



On the basis of our experience with spheres [171] and ellipsoids [31], we believe that our 
algorithm with rapid particle expansion produces final states that represent the MRJ state 
well. Here we use the largest possible initial 7 G (0.1 — 0.5) that is numerically feasible 
to ensure the generated superdisk and superball packings are maximally random jammed. 
We mainly focus on superdisks and superballs with deformation parameters p within the 
range 0.85 — 3.0, since extreme values of p associated with polyhedron-like shapes present 
numerical difficulties. 

All of the generated packings used in the subsequent analyses are verified to be at least 
collectively jammed using an "infinitesimal shrinkage" method 13], i.e., the particles in the 
packing are shrunk by a very small amount 39j and given random velocities. If no significant 
structural changes occur after the system "relaxes" after a sufficiently long enough time, the 
packing is considered to be collectively jammed. It is well established that the "infinitesimal 
shrinkage" method is robust, i.e., it always gives the same results for sphere packings as those 
obtained from a rigorous linear programming jamming test algorithm provided the amount 
of shrinkage from the jammed state is sufficiently small for a given number of particles within 



the periodic cell 



40]. 



7 





(a) p = 0.85 (b) p = 1.5 

FIG. 3: (color online). Typical configurations of MRJ packings of binary superdisks with different 
values of the deforn Le superdisks. 




FIG. 4: The density of MRJ packings of binary superdisks as a function of p. Insert: the average 
contact number Z as a function of p. 

A. Binary Mixtures of MRJ Superdisks 

In two dimensions, we study MRJ packings of a specific family of binary superdisk mix- 
tures in which the size ratio is k = 1.4 and the molar ratio is (3 = 1/3. The size ratio k is 
defined as the ratio of the diameter of large superdisks over that of the small superdisks; 
and the molar ratio /3 is defined as the number large superdisks over the number of small 
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(a) p = 0.85 



(b) p = 1.5 



FIG. 5: (color online). Typical configurations of MRJ packings of superballs with different values 
of the deformation parameter p. 

superdisks. We do not use monodispersed superdisk systems here because they are easily 
crystallized into ordered packings [25|. For p = 1, one obtains the binary circular-disk sys- 



tem which has been intensively studied as a prototypical glass former 41]. Typical jammed 



packing configurations are shown in Fig. [3j The density (ft and the average contact number 

lown in Fig. HJ which reveals that the initial rapid in- 



per particle Z function of p are s 

creases of (ft and Z are linear in \p — 1| [42j. The density (ft increases monotonically as p moves 
away from unity and shows no signs of a plateau, even for relatively large p. In addition, 
(ft quickly surpasses the density of the optimal binary circular-disk packing associated with 
the size and molar ratios employed here, which contains phase-separate regions of triangular 
lattice packings of different sized circular disks |41|. The quantity Z quickly reaches its 
plateau value Z* « 4.7 at p ~ 1.3, which is smaller than 2f = 6, indicating the packings are 
hypostatic. The cubatic order parameter P 4 is defined as P 4 = ((35 cos 4 9 — 30 cos 2 9 + 3)/8) , 
where 9 is the angle between the particle axis and the director, along which the principle 



axes of the particles have maximum mutual alignment 37]. The measured cubatic order 



parameter is P4 ~ 0.06 to 0.32 with the tendency to increase as \p — 1| grows. 
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FIG. 6: The density of MRJ packings of superballs as a function of p. Insert: the average contact 
number Z as a function of p. 



B. MRJ Packings of Monodisperse Superballs 

In three dimensions, mono dispersed superballs can be easily compressed into a jammed 
random packing due to geometrical frustration (i.e., the densest local particle arrangement 
cannot tile space). Typical jammed packing configurations are shown in Fig. [5j The quan- 
tities <p and Z as a function of p are shown in F ig. El respectively. As in two dimensions, 



4> and Z increase rapidly, in a cusp-like manner |43|, as the particles deviate from perfect 
sphere. increases monotonically as p moves away from unity, quickly goes beyond the 
optimal sphere packing density and shows no signs of plateau, even for relatively large p 
values. The contact number per particle Z reaches its plateau value Z* m 8.15 at p ~ 1.4, 
which is significantly smaller than 2f = 12, indicating that the packings are hypostatic. The 
measured cubatic order parameter is P4 ~ 0.03 to 0.21, which increases with \p — 1|. 



III. PACKING CHARACTERISTICS 



A. Rattlers 



MRJ packings generated in both two and three dimensions contain a small fraction of 
rattlers, i.e., particles that can wander freely within cages formed by their jammed non- 
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rattling neighbors. When p is close to unity, the fraction of rattlers is approximately 2.6% 
and 1.2% for two and three dimensions, respectively. As p moves away from unity, the 
fraction of rattlers decreases quickly and practically vanishes for large p (e.g., p > 2.75). 
This behavior results from the increasing protuberance of the particle shape, which makes it 
more difficult to form isotropic cages and also requires more average contacts per particle to 
achieve jamming. Note that rattlers are excluded when reporting average contact numbers 
in the following discussion. 



B. Packing Density 

The rapid increase of the density is mainly due to the broken rotational symmetry of 
the particles. In particular, the cubic-like (square-like) and octahedral-like particles are 
more efficient to cover the space than spheres (circular disks), i.e., near the jamming point 
the particles can rotate to accommodate the neighbors by orienting the "far corners" to 
fill the available gaps and thus cover more space. For small values of p, the increase in 
is also attributed to the expected increase in the number of contacting neighbors per 
particle, which means locally more particles can be packed in a given volume. The manner 
in which rotational symmetry is broken in superball packings is distinctly different from that 
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defined as the ratio of the radii 



in ellipsoid packings. For example, the asphericity 7 
of circumsphere and insphere of a nonspherical particle, is always bounded and close to unity 
for all values of p for superballs, while it can increase without limit as the largest aspect ratio 
a grows for ellipsoids. For very elongated or flake-like ellipsoids with large aspect ratios, 
the effect of a very anisotropic exclusion volume becomes dominant and causes the density 
of random ellipsoid packings to decrease. By contrast, the shape of superballs becomes 
more efficient in filling space as the deformation parameter deviates more from unity and 
thus results in a monotonically increasing density. The nonanalyticity of at p = 1 is 
also associated with the broken symmetry of superdisks and superballs. This nonanalytical 
behavior has also been observed in the optimal packings of these particles realized by various 



Bravais lattices 



25 



26]. This stands in contrast to the densest known ellipsoid packings, 



which are periodic packings with a two-particle basis possessing a smooth initial increase of 



as the aspect ratio moves away from unity 2l\ . 
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FIG. 7: (color online), (a) Distribution of contact numbers for different p values for MRJ packings 
of superdisks (upper panel) and superballs (lower panel), (b) Local packing structures with more 
contacts than average (shown in blue) and those with less contacts than average (shown in pink) 
in two-dimensional superdisk packings for different p values. 

C. Hypostaticity and Nongeneric Local Structures 



There have been conjectures 15|, ll6] that frictionless random packings have just enough 
constraints to completely statically define the system (i.e., it is isostatic), i.e., for large 
packings, one has Z = 2f. It has been shown both experimentally and computationally 



that although the isostatic conjecture [15( holds for large sphere packings 



40 



31 



44j, it is 



32 1 . It was found 



generally not applicable to nonspherical particles, such as ellipsoids 
that even for ellipsoids with large aspect ratios, Z is still slightly below 2f 31]. 

Here we observe that in MRJ packings of superdisks and superballs Z is significantly 
smaller than 2f for all values of p examined, i.e., the packings are significantly hypostatic. 
The hypostatic packings result from the competition between fx translational and ro- 
tational degrees of freedom of the particles (/ = //? + fr) m developing the contacting 
networks close to the jamming point. In particular, although it is true that to constrain 
the translational degrees of freedom each particle needs at least 2Jt contacts, rotational 
degrees of freedom can be blocked with less than 2ff> additional contacts per particle if 
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the curvatures at the contacting points are sufficiently small [32|. In addition, due to the 
relatively small asphericity 7 of superdisks and superballs, there is little reason to expect the 
rotational motions of these particles (especially those with p close to unity) would be frozen 
even when they are translational trapped and may only rattle inside small "cages" formed 
ibors. Near the jamming point, it is expected that the particles can rotate 



by their neig 

significantly 45j until the actual jamming point is reached, at which rotational jamming 
will also come into play, and rotational degrees of freedom are frozen with the number of 
additional contacts much less than 2/r. This is in contrast to hypostatic MRJ packings 
of ellipsoids with large aspect ratios, for which the translational and rotational degrees of 
freedom are on the same footing and, thus, the average contact number per particle is only 
slightly below twice the number of total degrees of freedom. 

Furthermore, the local geometry of the MRJ packings is necessarily nontrivially correlated 
(nongeneric), i.e., all the normal vectors at the points of contact for a particle should intersect 
at a common point to achieve torque balance and block rotations. In light of the isostatic 
conjecture, the local packing structures are less nongeneric when they possess larger contact 
numbers so that the constraining neighbors are less correlated. The truly generic local 
packing structures should have Z = 2f per particle, for which the constraining neighbors 
could be completely uncorrelated. To characterize the "nongenericity" of the packings, we 
compute G ng , the fraction of local structures composed of particles with less contacts Zi oca i 
than average Z average , i.e., 

GN(Z local ^ Z avera g e ^j 
ng = • \ Z ) 

^ total 

A larger G ng indicates a larger degree of nongenericity. We find G ng is approximately 
0.65 in two dimensions and 0.78 in three dimensions when p is close to unity, which quickly 
decreases and plateaus at 0.6 and 0.68, respectively as p increases. Figure [7] shows the 
distribution of contact numbers for different p values and the topology of the local structures 
contributing to G ng . It can be seen that as p moves away from unity, the distributions become 
more skewed as the means shift to larger Z. Moreover, the subset of particles associated with 
the nongeneric structures do not percolate. We do not observe any tendency of increasing 
Z even for the largest p values that are computationally feasible and we expect that MRJ 
packings in the cubic limit are also hypostatic. It is noteworthy that isostatic random 
packings of superdisks and superballs are difficult to construct, since achieving isostaticity 
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(a) (b) (c) 

FIG. 8: (color online). Nongeneric locally jammed configurations associated with three fixed su- 
perdisks (pink) and the trapped one (blue). In each configuration, the central superdisk is approx- 
imately aligned with one of its fixed neighbors to form contacts associated with small curvatures 
to block rotation. 



requires Z — 2f = 12 which is necessarily associated with translational crystallization [461 ] . 

We note that the aforementioned nongeneric structures (see Fig. [3(b)) are not rare. In 
particular, a nonspherical particle can be rotationally jammed if it has neighbors that can 
translationally jam the particle [32)]. To illustrate this point, we will consider a small packing 
composed of four superdisks in two dimensions. Now we show that one can locally jam a 
superdisk by three contacting neighbor superdisks. Translational jamming requires that 
the centroids of the neighbors cannot lie in the same semi-circle around the centroid of 
the central superdisk. Suppose a superdisk is translationally trapped (not jammed) by its 
three neighbors, whose positions and orientations are fixed. This four-particle configuration 
has four degrees of freedom: two translational and one rotational degrees of freedom of the 
trapped particle as well as the expansion of the particles. To obtain a jammed configuration, 
the four degrees of freedom need to be completely constrained. This can be achieved by the 
three contact conditions for the jammed particles and its neighbors and the requirement 
that the three inward normal vectors at the contacting points meet at a common point, a 
sufficient condition for torque balance 



321 ] . Thus, one has four independent equations for 
the four degrees of freedom; see the Appendix for details. 

Figure M shows the nongeneric jammed configurations associated with three specific fixed 
trapping superdisks. The multiplicity of the configurations is due to the multiple solutions of 
the equations. The jamming configurations can be also obtained using the DTS algorithm, 
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which allows the trapped particle to translate and rotate and allows all the particles to grow. 
Although the above analysis is for local jamming (i.e., the neighbors of the central particle 
are fixed), it is reasonable to expect that collective particle rearrangements further facilitate 
the formation of nontrivial orientational correlations and thus, enable a larger number of 
nongeneric jamming configurations. Indeed, the numerous hypostatic jammed packings that 
we found from our simulations strengthen our argument that nongeneric structures are not 
rare. 



D. Nonvanishing Orientational Order 

We also observe the increase of the orientational order (measured by P4) associated with 
the increasing p values, although the largest possible expansion rate 7 has been used to 



suppress the formation of orders (i.e., to maintain the maximal degree of randomness) [47]. 
As p deviates from unity, the particle shape develops "edges" and "corners" with large 
curvatures, which may not be able to block rotational unjamming motions if contacts occur 
at significantly curved regions of the particle surface. On the other hand, low-curvature 
contacts are more favorable, which is associated with partial alignments of the particles. The 
tendency of particle alignments to form low-curvature surface contacts required by jamming 
becomes stronger as the particle moves further away from the sphere point. Thus, there is 
also a competition between orientational disorder and jamming for packings of superdisks 
and superballs, resulting from their unique symmetry-breaking manner, which has not been 
observed in random packings of ellipsoids. Due to numerical difficulties, we could not use the 
DTS algorithm to study the random jammed packings of particles with extreme shapes, i.e., 
in the limit p — > 0.5 and p — > 00. However, it reasonable to expect considerable orientational 
ordering in such packings. 



IV. CONCLUSIONS 



In this paper, we studied the maximally random jammed packings of superdisks and su- 
perballs. The packing densities increase dramatically and nonanalytically as one moves away 
from the circular-disk and sphere point (p = 1) and the packings are hypostatic. To achieve 
jamming, the local arrangements of particles are necessarily non-trivially correlated and we 
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term these structures "nongeneric" in light of the correlations. The degree of "nongenericity" 
of the packings is quantitatively characterized by the fraction of local structures composed 
of particles with less contacts than average. Moreover, we showed that such seemingly "spe- 
cial" packing configurations are not rare. As the anisotropy of the particles increases, the 
fraction of rattlers decreases while the minimal orientational order increases. The novel fea- 
tures arising in MRJ packings of superdisks and superballs result from the unique manner 
in which rotational symmetry is broken. This makes such packings distinctly different from 
other known MRJ packings of nonspherical particles, such ellipsoids and ellipses. 

The ability to produce dense random packings using superballs casts new lights on sev- 
eral industrial processes, such as sintering and ceramic formation, where interest exists in 
increasing the density of powder particles to be fused. If superball-like particles instead 
of spherical particles are used, the packing density of a randomly poured and compacted 
powder could be increased to a value surpassing that of the maximal sphere-packing density. 
We note that superdisks and superballs can be experimentally mass produced using current 
lithography techniques. Understanding the statistical thermodynamics of the jamming tran- 
sition of superdisks and superballs, especially the role of rotational and translational degrees 
of freedom for different deformation parameters is a subject that merits future investiga- 
tion. Such studies could may our understanding of the nature of glass transitions, since the 
preponderance of previous investigations have focused on spherical particles. 



APPENDIX: EQUATIONS FOR LOCALLY JAMMED FOUR-SUPERDISK 
CONFIGURATIONS 

In this section, we provide the equations that determine locally jammed four-superdisk 
configurations composed of a trapped central particle and three fixed contacting neighbors 
(see Fig. |HJ). In particular, the boundary of a superdisk with radius R is define by 



R 



2p 



+ 



X2 

R 



2p 



(A-l) 



which can be also expressed by the parametric equations 



x\{9) = | cos 9\ p ■ R ■ sign(cos 9), 
#2(0) = I sin#|p ■ R ■ sign(sm9), 



(A-2) 
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where sign(x) gives the sign of argument x. Let the centroids of the three fixed neighbors 
be (<Zj, bi) (i = 1, 2, 3), and the orientations be 9 { . Their boundaries are then given by 



x[ (9) = Of + cos^l cos#|p • R ■ sign(cos9) + sin 6*;| sin#|p • R ■ sign(sm9), 
x 2 = bi — sin ^ | cos#|p ■ R ■ sign(cosd) + cos^l sin#|p • R ■ sign(sin9). 

Similarly, if the centroid of the central particle is at (a Q , b ) and its orientation is characterized 
by 9 Q , its boundary is specified by 



x°(9) = a Q + cos 9 Q \ cos 9\ p ■ R ■ sign(cos9) + sin# | sin^lp ■ R ■ sign(sm9), 
x%(9) = b Q — sin^ | cos 6\p ■ R ■ sign(cos9) + cos^ G | sin9\p ■ R ■ sign(sm9). 



(A-4) 



Since the positions and orientations of the three neighbors are fixed, the four-particle system 
has four degrees of freedom, namely the position (a Q , b Q ) and orientation 9 Q of the central 
particle, as well as the radius R of all particles, as discussed in Sec. III.C. 

In the jammed configuration, the central particle contacts all its three neighbors. From 
Eqs. f lA-ll) and (1A-41) . the contact point (xil,x^) between neighbor particle % and the central 
particle can be expressed in terms of (a ,b ,9 , R), which must also lie on the boundary of 
the neighboring particle i, i.e., 



.(0. 

'lc 1 



K, &o, R) — a i 



R 



2p 



a>o, b , 9 a , R) — hi 



R 



2p 



(A-5) 



for i = 1,2,3. This leads to three equations in the variables a ,b ,9 , and R. In addition, 
to achieve jamming the three normals at contacts must meet at a common point, which 
guarantees torque balance. The normals at contacts are along the lines 



(x 2 - 4c) = - 



dxf/d9 



dxf/dB 



{xi-xf c ). 



(A-6) 



The aforementioned torque balance condition requires that the three lines given by Eq. (|A-6[) 
must intersect at a common point. This leads to another equation in the variables a D , b a , 9 Q , 
and R. Thus, there are four independent equations for the four degrees of freedom and 
(a , b , 9 Q , R) can be completely determined for a locally jammed four-particle configuration. 
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